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ABSTRACT 


This  report  contains  the  results  of  an  analytic  investigation  of  the 
Transverse  Doppler  Pattern  Measurement  Technique  for  determining 
the  far  field  radiation  characteristics  of  large  antennas  from  near  field 
measurements.  The  technique  is  based  on  using  Doppler  signal  process¬ 
ing  to  obtain  the  information  necessary  for  predicting  the  far  field  pat¬ 
terns  of  large  antennas  in  their  site  environments.  Primary  emphasis 
has  been  placed  on  establishing  a  rigorous  mathematical  description  of 
the  measurement  process  from  which  system  parameters  and  perform¬ 
ance  may  be  determined. 

Several  mathematical  models  have  been  developed  which  differ 
both  in  the  analytical  approach  and  in  the  physical  make-up  of  the  meas¬ 
uring  system.  Mathematical  approximations  for  the  diffraction  field  of 
the  aperture  have  been  expanded  in  several  coordinate  frames.  The  sig¬ 
nal  processing  systems  investigated  include  multichannel  processing, 
synthetic  aperture  processing,  and  directional  and  nondirectional  sam¬ 
pling  antennas.  Attention  has  been  given  to  linear  arrays,  rectangular 
arrays  with  separable  and  nonsepar  able  distributions,  circular  apertures 
with  circularly  symmetric  distribution,  and  arbitrarily  shaped  radiating 
apertures  with  large  linear  phase  deviations.  Two  experimental  proce¬ 
dures  have  been  outlined  for  the  purpose  of  verifying  the  feasibility  of 
the  analytical  work. 
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I. 


INTRODUCTION 


The  far  field  behavior  of  antennas  is  customarily  measured  at  dis- 

2 

tances  which  exceed  or  are  at  least  equal  to  2D  /V,  where  D  is  the 
aperture  diameter  and  X.  is  the  operating  wavelength  in  free  space.  For 
many  large  antennas,  especially  those  used  for  radio  astronomy,  con¬ 
ventional  methods  of  pattern  measurement  at  this  distance  are  always 
difficult  and  often  impossible.  For  this  reason,  various  other  techniques 
have  been  utilized  in  obtaining  desired  data  on  the  far  field  performance 
of  these  large  antennas.  These  techniques  include  actual  far  field  meas¬ 
urements  using  aircraft,  balloons,  and  radio  stars  to  obtain  the  required 
measurement  distance  as  well  as  near  field  measurements  such  as  prob¬ 
ing  the  aperture  field  to  obtain  the  aperture  distribution  and  focusing  the 

12  3 

antenna  within  the  near  field.  ’  ’  Such  techniques  are,  in  general, 

straightforward  extensions  of  conventional  pattern  measurement  proce¬ 
dures  and  are  more  or  less  satisfactory  depending  on  the  particular 
antenna  being  evaluated.  They  are  usually  involved  and  time  consuming 
and,  in  some  cases,  require  modification  of  the  antenna  from  its  opera¬ 
tional  configuration  with  a  resultant  uncertainty  in  the  measured  data. 

Described  in  this  report  is  the  analytic  investigation  of  the  appli¬ 
cation  of  Doppler  techniques  to  the  large  antenna  pattern  measurement 
problem  in  an  attempt  to  overcome  the  various  difficulties  associated 
with  more  conventional  methods  of  measurement.  The  aperture  distri¬ 
bution  is  obtained  by  processing  the  signal  received  by  a  probe  antenna 
which  moves  at  high  velocity  across  the  radiating  aperture  at  a  distance 
which  is  a  small  fraction  of  its  far  field  distance.  The  far  field  radia¬ 
tion  pattern  is  then  obtained  directly  from  the  measured  aperture  distri¬ 
bution  by  performing  the  radiation  integration.  The  analytic  expressions 
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which  form  the  basis  of  the  Transverse  Doppler  Pattern  Measurement 
Technique  are  considered  in  detail  and  various  methods  of  signal  proc¬ 
essing  are  examined.  The  technique  is  applied  to  various  antenna  con¬ 
figurations  and  computations  of  the  results  are  presented.  Experiments 
for  determining  the  feasibility  of  the  measurement  technique  are  also 
discussed. 
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II-  GENERAL  AKIAL^SIS 


Under  the  ass.  umption  that  the  antenna  whose  characteristics  are 

to  be  measured  may  be  xepresented  by  an  aperture  lying  in  the  plane 

and  that  the  apertur-  e  field  is  uniformly  polarized  in  the  ^-direction,  the 

4 

diffraction  field  is  -ggiven  by  (using  Silver's  scalar  field  approximation  ): 


U 


exp  (jWQt) 
4rr 


£ 


exp  (-jkr) 


(ik+;)  VJ1  +  jkVS_ 


d^dr| 


for  the  geometry  of  Figure  la.  In  the  equation, 

F(£,q)  =  A(|  »  rj)  exp  Q'i'd,  q  )J,the  aperture  distribution 

$  is  a.  unit  vector  in  the  direction  from  the  aperture 
poiimt  4,  r]  to  the  field  point  P 

§  is  a*,  unit  vector  along  a  ray  through  the  aperture. 

Equation  (1)  may  b^s  considerably  simplified  under  the  following  condi¬ 
tions: 


k  >>  —  or  kr  >>  1 
r 


A  A  A  A 

u  _  i»r  ai»R  =  cos  0  a  1 
z  1  z  1 


iii  __  i  s  1 
z 


The  first  condition  requires  that  the  measurement  distance  be  large  com¬ 
pared  with  the  operating  wavelength.  The  second  requires  that  the 
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measurement  distance  be  large  enough  so  that  the  unit  vector  r  , 
directed  from  an  arbitrarily  located  current  element  on  the  aperture  to 
the  field  point  P  ,  is  essentially  parallel  to  the  unit  vector  &  directed 
from  the  center  of  the  aperture  to  the  field  point,  and  that  the  antenna 
main  beam  be  at  or  near  broadside.  The  third  condition  requires  that 
the  phase  distribution  be  essentially  planar.  These  conditions  are,  in 
general,  adequately  met  for  large  antennas  of  the  type  being  considered 
here,  and  the  scalar  diffraction  equation  reduces  to 


U 

P 


jkexp  (j«Qt) 
2tr 


m.r,) 


exp  ( -jkr ) 
r 


d£dr| 


(2) 


Now  assume  that  the  diffraction  field  is  sampled  with  an  omnidirectional 
probe  antenna  which  is  moving  across  the  aperture  at  a  constant  height 
h  with  a  constant  velocity  v  .  The  probe  path  is  described  by  the  space 
coordinates  x  =  0,  y  is  variable,  and  z  =  h.  The  y-axis  is  coincident 
with  the  q-axis,  and  the  origin  of  the  space  coordinate  system  coincides 
with  that  of  the  aperture  coordinate  system  as  indicated  in  Figure  lb. 


In  rectangular  coordinates,  r  may  be  expressed  as 


r 


+  (  t)  -  y) 


2 


1/2 


(3) 


Thus,  in  terms  of  the  parameters,  h  and  y,  the  radiated  field  along 
the  probe  path  is 
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dem.ro 


gy(h»y)  = 


jk  exp  (ju> 


2tt 


C2(ti) 


c^r)) 


exp 


(-jk[h 


u2  ,2  .  .2 

h  +  £  +  (y  -  r| ) 


D 


[h2  +  i2  +  (y  -  T))2] 


1/2 


For  y  a  linear  function  of  time  (corresponding  to  a  constant  probe 
velocity) 


y  =  vt  ;  -oo<t<+oo 


and 


jk  exp  (jw  t) 

,  gy(h,vt)  =  - —2. 


0^2  PC2(Tl) 

'  ^  \(r\) 


d|F(g,r,) 


exp{.Jk[,vt-„2th2t|2]‘/2} 


[(vt-,)2+h2  +  e2] 


1/2 


Now  let  us  consider  the  effect  of  passing  this  time  dependent  signal 
through  a  bandpass  filter  as  it  is  received  at  the  probe  antenna.  The 
spectral  density  of  the  signal  after  passage  through  the  filter  is  found 
using  the  Fourier  transform  pair: 


1 


g(t) 


G(u) 


n|  2tt 


1 

^  2tr 


G(w)  exp  (jwt)  dw 


g(t)  exp  (-jwt)  dt 


(6) 


For  a  filter  with  a  characteristic  function  H(w)  ,  the  spectral  density  of 
the  output  signal  is 


S  (w)  =  H(w)  G(w) 

and  the  output  signal  is  (using  Equation  (6)  above) 


(7) 


S  (t) 


i  C00 

nJ  2tt  ^-00 


S  (w)  exp  (jwt)  dw 


p2  rczM  p« 

\  dq  \  d|F(|,ri)  \  * 

Uti1  Hr^q)  J-oo 


(2ir)  ^ 


dwH(w)  exp  (jwt) 


ioo 

•oo 


f-jk|(vt*  -  T,)2  +h2  +  £2]  } 

- m -  exp  [-j(w  -  wq)  t]  dt' 


,'oo  exp  j ~jk|(vt'  -  tj) 
•  \  k 


fvt'  -  t,)2  +h2  +  e2] 


(8) 


This  equation  is,  in  general,  quite  difficult  to  evaluate  because  of  the 
complicated  way  in  which  both  the  space  and  time  variables  appear. 
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Some  simplification  of  the  mechanics  of  the  integrals  may  be  obtained 
by  carefully  choosing  the  filter  characteristics 


As  examples,  consider  the  filter  functions 


r 


=  ■{ 


B  B*^ 

1  ;  w  -  —  <  oj  <  oj  +  — 

o  2  o  2 


0  ;  elsewhere 


(9 


„  ,  .  SinV  2B 
H  (w)  =  — 

L  CO  -  CO 


CO  -  CO 


(II 


where  B  is  the  filter  bandwidth.  For  these  cases  Equation  (8)  may  be 
integrated  to  give 


jk exp  (jw  t)  0^2  PC2^ 

V‘>  ■  -  2  -  J  <"■  J  ,  . 

1  2  IT  C^Tl) 


d|F(e,T!) 


i'oo 

- 


exp 


{[v(t"  U|  -  r|]2  +  h2  +  i2} 


“"f  *" 

t" 


dt" 


for  the  flat  top  bandpass  filter  and 
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jkexp(jw  t)  0^2  PC2^ 

S  (t)  = - - - 2-  \  dT1  \  d€F(£,q) 

V2  4,r  Jr,1  JCi(t!) 


for  the  Sln  U  filter.  However,  it  appears  that  neither  of  these  equations 
can  be  integrated  further  in  closed  form. 

At  this  point  it  is  advisable  to  consider  the  physical  counterpart  of 
the  mathematical  expression  derived  above.  For  a  vehicle  moving  at  a 
velocity  v  and  a  height  h  above  a  plane  surface  which  is  uniformly  radi¬ 
ating  energy  at  a  frequency  f^  ,  the  frequency  of  the  energy  received  at 
the  vehicle  from  any  point  on  the  plane  will  be  shifted  from  f^  by  an 
amount  proportional  to  the  relative  velocity  with  which  the  vehicle  is  ap¬ 
proaching  or  receding  from  the  radiating  point.  This  is  the  well-known 
Doppler  frequency  relation.  It  may  be  readily  shown  that  the  points  on 
the  plane  for  which  the  received  energy  has  a  constant  Doppler  shift  at 
any  given  instant  form  a  family  of  hyperbolas  with  a  straight  line  directly 
below  the  vehicle  and  transverse  to  its  path  as  the  limiting  case  corre¬ 
sponding  to  zero  Doppler  shift.  Therefore,  filtering  at  fQ  allows  passage 
of  only  those  frequency  components  which  correspond  to  signals  radiated 
from  the  region  between  two  hyperbolas  on  the  plane  surface  as  indicated 
in  Figure  2a.  This  sampling  region  moves  with  the  vehicle  as  it  passes 
-over  the  plane.  If  the  portion  of  the  plane  which  radiates  is  now  con¬ 
strained  to  be  an  antenna  aperture,  the  effect  is  that  of  moving  a  samp¬ 
ling  function  progressively  over  the  aperture  as  shown  in  Figure  2b. 
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Figure  2b.  Doppler  Sampling  Region  on 
a  Planar  Aperture 


Using  this  physical  picture,  a  simplification  of  Equation  (12)  can 
be  made  by  recognizing  that  the  limits  of  integration  in  t"  are  related 
to  the  filter  time  constant  and  that  the  time  constant  is  proportional  to 
the  sampling  region  width  divided  by  the  probe  velocity.  For  practical 
values  of  the  sampling  region  width  (1/10  to  1/100  of  the  aperture  dimen¬ 
sions),  the  integrand  can  be  approximated  in  a  manner  which  is  similar 

gyn  |  «ilr  j*  | 

to  the  customary  far  field  approximation  of  — ■  —  - - ,  In  the  phase 

factor,  let 


jjv(t»  +  t)  -  n]2  +  h2  +  52}'/a  »  (h2  +  e2)'/2 


1  + 


vtn(vt  -  11  ) 
(h2  +  iz) 


(1 


and,  in  the  amplitude  factor,  let 


-a/2 


("2  *  e2) 


aO/2 


{[v(f  + 1)  -  q]2  +  h2  +  e2 
With  these  substitutions,  the  integration  in  t"  can  be  performed  to  give 


(1 


S  (t)  = 


=  zV  exp(jV)  J 


"z  r'a'11' 

dq  \ 

T1L  ^c^q) 


d|FU,T,) 


sin 


f  2  2V^2 

•  exp  -jk\Ji  +  £  J 


~kv 

(vt  - 

n)  1 

2B  / 

r 

,v/2 

( 

,h 

7  J 

kv(vt  -  T|) 


(1 
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The  result  is  consistent  with  the  physical  operation  of  the  system.  The 
-■1^  U  function  in  the  integrand  is  sharply  peaked  at  vt  =  t)  and  has  the 
effect  of  moving  a  sampling  function  across  the  aperture  with  a  velocity 
v.  The  resolution  of  the  sampling  function  is  dependent  on  the  param¬ 
eters  k,  v,  B,  and  h.  The  transverse  aperture  variable  £  also  appears 
both  in  the  sampling  function  and  in  the  phase  factor.  In  the  sampling 
function,  the  effect  is  an  aperture  resolution  which  is  a  function  of  the 
transverse  variable.  This  results  from  the  fact  that  the  lines  on  the 
aperture  plane  corresponding  to  constant  Doppler  frequency  are  hyper¬ 
bolas  rather  than  straight  lines.  The  appearance  of  the  variable  £  in 
the  phase  factor  is  caused  by  the  fact  that  the  probe  antenna  is  in  the 
near  field  of  the  transverse  dimension  of  the  sampling  strip.  This  de¬ 
pendence  is  not  required  in  the  q  dimension  since  the  probe  is  in  the 
far  field  for  the  width  of  the  sampling  strip  in  q  . 

To  summarize,  within  the  accuracy  implicit  in  the  approximations 
of  Equations  (13)  and  (14),  Equation  (15)  represents  the  output  signal 
from  a  rapidly  moving  Doppler  antenna  system  as  it  passes  over  a  gen¬ 
eral  two-dimensional  planar  aperture.  It  has  been  assumed  that  the 
Doppler  filter  has  S1^'U  bandpass  characteristics,  and  that  the  probe 
antenna  has  an  omnidirectional  lower  hemisphere  coverage  and  is  -lov¬ 
ing  with  a  constant  velocity  v  at  a  constant  height  h  . 
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III.  SPECIFIC  APPLICATIONS 


A  large  number  of  practical  antennas  may  be  represented  mathe¬ 
matically  by  three  planar  configurations:  (1)  the  line  source,  (2)  the 
rectangular  aperture  with  its  aperture  distribution  separable  in  the 
orthogonal  planes  containing  the  major  dimensions  of  the  aperture,  and 
(3)  the  circular  aperture  with  a  circularly  symmetric  aperture  distribu¬ 
tion.  The  Transverse  Doppler  Pattern  Measuring  Technique  has  been 
applied  to  these  configurations,  and  the  results  of  the  analytic  investiga¬ 
tions  are  presented  below. 

A.  Line  Source  with  Uniform  Phase 


Perhaps  the  simplest  aperture  to  consider  is  the  line  source 
with  uniform  phase.  In  addition,  a  detailed  study  of  the  application  of  the 
measurement  technique  to  the  line  source  provides  valuable  information 
about  basic  system  parameters. 


to 


For  a  line  source  of  length  a  in  q  ,  Equation  (15)  reduces 


S  (t)  =  exp  (jwQt) 


dq  A(  q) 


kv 

2Bh 


(vt  -  q) 


kv(vt  -  q) 


(1 


where  a  constant  phase  factor  and  a  multiplicative  constant  have  been 
omitted.  The  aperture  sampling  function  is 

sin 

kv(vt  -  q) 


kv 

2Bh 


(vt  -  q) 
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Determination  of  the  optimum  values  of  the  parameter  — requires 
that  Equation  (16)  be  evaluated  for  the  particular  aperture  distribution 
which  is  being  measured.  However,  intuitive  considerations  can  be  used 
to  approximate  the  parameters.  If  the  aperture  dimension  is  a,  the 
resolution  on  the  aperture  can  be  expressed  as  a/N,  where  N  is  the 
number  of  sampling  widths  within  the  aperture.  If  the  sampling  resolu¬ 
tion  is  chosen  to  be  the  3  db  width  of  the  sampling  function,  the  following 
relation  of  system  parameters  is  obtained. 


kva 

5.6  BhN 


(18 


An  additional  constraint  on  the  system  parameters  is  im¬ 
posed  by  the  approximations  made  in  the  evaluation  of  the  t"  integral 
in  Equation  (12).  The  "far  zone"  approximation  of  the  integrand  re¬ 
quires  that 


(19 


or 


h  = 


a  >  1 


(20 


When  this  condition  is  substituted  in  Equation  (18),  the  equation  takes 
the  form 


v 

B 


(21 
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Therefore,  the  two  equations  governing  the  choice  of  the  system  param¬ 
eters  are 


kh  = 


(2 


kv  _  ka  r  7T  [ 
B  =  N  V.  5.6  a  J 


(2 


where  the  form  of  Equations  (20)  and  (21)  have  been  modified  for  conven¬ 
ience  in  graphing.  These  equations  are  plotted  in  Figures  3  and  4  with 
a  as  a  parameter.  The  use  of  these  curves  to  determine  the  operational 
parameters  of  the  measuring  system  can  be  illustrated  by  a  typical 
sample  problem.  Consider  a  line  source  1000  feet  in  length  operating 
at  a  wavelength  of  2  feet: 

ka  =  3142  and,  for  a  linear  resolution  corresponding  to 
N  =  10,  ka/N  =  314.2.  From  Figure  3  with  a  =  2,  kv/B  =  88.3.  Thus, 
v/B  =  28.1.  If  v  =  600  mph  =  880  ft/sec,  B  =  31  cps.  From  Figure  4, 
for  kv/B  =  88.3  and  a  =  2,  kh  =  4960,  and  h  =  1580  feet. 


Summarizing, 

h  =  1580  feet 
v  =  880  ft/sec 
B  =  31  cps 
N  =  10 


It  should  be  noted  that  the  same  value  of  a  must  be  used  for  both  graphs. 
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1000 


ka 

N 


Figure  3. 


ic  v  kci 

Relation  Between  — —  and  —  for 
B  N 


a  >  1 


16 


17 


The  above  example  indicates  that  the  system  designer  has  a 
considerable  amount  of  latitude  in  choosing  the  system  parameters  pro¬ 
vided  the  parameters  are  selected  consistent  with  Equations  (22)  and 
(23). 


To  determine  what  effect  the  choice  of  the  composite  param- 
kva 

eter  K  =  ■  has  upon  the  accuracy  with  which  the  aperture  distribution 

2Bh 

is  reproduced.  Equation  (16)  has  been  evaluated  numerically  in  the  form 


S(T) 


Pl/2 
\  A(V) 
'-'-1/2 


sin  [k(t  -  y)1 
K(t  -y) 


dv 


(24 


for  various  typical  uniform  phase  line  source  distributions  and  values 

of  the  parameter  K.  The  results  of  these  calculations  are  presented 

5 

in  Scientific  Report  No.  1  on  this  contract  on  pages  19  through  23.  It 
can  be  seen  that  the  accuracy  with  which  the  aperture  distribution  is 
measured  increases  as  the  value  of  K  increases.  This  is  as  expected 
since  increasing  the  value  of  K  has  the  effect  of  narrowing  the  sampling 
function  and  correspondingly,  the  sampling  width.  In  fact,  it  may  be 
easily  shown  that,  if  the  sampling  function  is  the  Dirac  delta  function 
(the  case  for  K — ►  »)  ,  the  aperture  distribution  is  measured  exactly. 

It  is  also  apparent  that  the  best  reproduction  is  obtained  for  those  dis¬ 
tributions  which  are  zero  at  the  ends  of  the  aperture  rather  than  for  those 
which  have  abrupt  discontinuities  at  the  ends.  This  is  to  be  expected 
from  the  analogy  with  Fourier  representations  of  abruptly  discontinuous 
functions. 
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B.  Rectangular  Apertures  with  Separable 
Distributions  and  Uniform  Phase 


If  the  conditions  for  a  rectangular  aperture  with  a  uniformly 
phased  separable  distribution,  i.  e.  , 

F(£,t,)  =  f(€)  f(r,> 

Cjft,)  =  -b/2 

c2(ii)  =  b/2 
r]l  =  -a/2 

^2  =  a/2 

are  introduced  in  Equation  (15),  it  becomes 


kv(vt  -  ti) 


It  should  be  noted,  as  mentioned  above,  that  the  transverse  variable  £ 
appears  in  both  the  phase  factor  and  the  sampling  function.  However, 
for  typical  values  of  h  and  £,  the  variation  in  the  sampling  function 
resolution  is  less  than  one  percent  over  the  full  excursion  in  £  and  may 
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be  neglected  with  very  little  error.  With  this  approximation  Equation 
(25)  may  be  written 


Sy(t) 


Jfc 

2tt 


exp  (jwQt) 


Ob/2 

\  exp 

vJ-b/2 


dqf(q) 


kv 

sin 

2Bh  (Vt  " 

kv(vt  -  q ) 


(2< 


Thus  the  dependence  in  5  appears  only  as  a  constant  which  modifies  the 
amplitude  of  the  expression.  Therefore  the  discussion  in  the  preceding 
section  with  regard  to  linear  arrays  may  be  applied  directly  to  the  con¬ 
sideration  of  rectangular  apertures  with  separable  distributions.  A 
second  measurement  orthogonally  oriented  in  space  (with  the  roles  of 
£  and  q  in  Equation  (25)  interchanged)  is  required  to  obtain  the  distri¬ 
bution  in  %  . 


C.  Nonuniformly  Phased  Linear  and  Separable, 

Rectangular  Distributions 

If,  in  Equation  (16),  the  aperture  excitation  function  is  as¬ 
sumed  to  be  complex,  the  equation  has  the  form 


Sy(t)  =  exp  (jwQt) 


Pa/2 

J  d  ti  A(  T| )  exp  jl  <j>  ( q  jj 


sin 


kv 

2Bh 


(vt  -  q) 


kv(vt  -  q) 


(27 


where  F(q)  =  A(q)  exp jj <j> ( q jj .  (A(q)  and  <f>(q)  are  real.)  In  partic¬ 

ular  for  the  usual  linear  phase  shift  required  for  scanning, 
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4>(*i)  =  P*i 


(2 


with  p  a  constant.  This  equation  has  been  evaluated  numerically  in  the 
form 


ru 

S(t)  =  \ 

^-1/ 


2  sin 

A(y)  exp  (jpY) 

2 


K(t  -  y)~ 
K(t  -  y) 


dy 


(2 


for  two  aperture  distributions,  uniform  and  (1  -  2y  )  ,  and  values  of  p 

corresponding  to  a  total  phase  variation  across  the  aperture  of  1,  5,  10, 

and  20  radians.  The  results  of  the  computations  are  presented  on  pages 

28  through  36  of  Scientific  Report  No.  1  on  this  contract.^  In  general, 

the  reproductions  of  the  amplitude  compare  favorably  with  those  obtained 

for  the  in-phase  apertures,  and  reproduction  of  the  phase  is  excellent 

over  the  aperture.  The  exceptions  occur  for  small  values  of  K  and 

P  =  20.  The  relation  of  these  parameters  to  the  physir  system  may  be 

kva 


determined  from  a  consideration  of  Equation  (18).  Since  K  = 
have 


2Bh  ’ 


we 


kva  K 

5.6  BhN  =  1  =  2.8  N 


or 


N  = 


K 

2.8 


(3 


where  N  is  the  number  of  sampling  widths  within  the  aperture.  There¬ 
fore,  a  value  of  K  =  10  corresponds  to  a  sampling  strip  approximately 
0.28  of  the  aperture  dimension  in  width.  Since  the  aperture  phase  slope 
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is  20  radians,  there  is  a  linear  phase  shift  of  5.6  radians  or  almost  a 
full  cycle  across  each  sampling  width.  It  is  not  surprising,  therefore, 
that  the  phase  and  amplitude  resolution  is  poor  for  these  values  of  the 
parameters. 


The  above  values  of  phase  variation  over  a  large  aperture 
correspond  to  relatively  small  scanning  of  the  antenna  beam.  For  ex¬ 
ample,  for  a  500X  aperture,  (3  =  20  results  in  a  scan  of  approximately 
two  beamwidths.  For  scanning  to  an  angle  of  20°  (of  the  order  of  120 
beamwidths),  a  value  of  p  a  1075  is  required.  Numerical  evaluations 
of  Equation  (31)  for  this  value  of  (3  and  several  distribution  parameters 
were  also  performed.  These  are  definitely  unsatisfactory  with  regard 
to  both  phase  and  amplitude.  As  discussed  above,  however,  this  is  not 
unexpected.  Even  for  K  =  60  ,  corresponding  to  approximately  21.4 
sampling  widths  within  the  aperture,  there  is  a  phase  slope  of  approxi¬ 
mately  50  radians  across  a  sampling  strip.  Since  the  signal  out  of  the 
filter  is  essentially  an  average  signal  from  the  sampling  strip,  it  can  be 
shown  that  the  amplitude  of  the  signal  tends  to  be  zero  except  near  the 
edges  of  the  aperture  where  a  nonsymmetrical  condition  exists  with  re¬ 
spect  to  the  sampling  function. 


22 


1 


D.  Circular  Aperture  with  a  Circularly  Symmetric 
Aperture  Distribution 


For  a  circular  aperture  of  radius  a/ 2  with  a  circularly 
symmetric  aperture  distribution,  Equation  (15)  has  the  form 


•k  Ca/2 

S  (t)  =  p  exp  (j«  t)  \  dq 

V  2,r  °  J-a/2 


deF(e.ri) 


sin 

kv  (vt  -  q) 

_2B  (h2*e2r. 

kv(vt  -  q  ) 

(32 


where  F(§,q)  is  circularly  symmetric  in  %  and  q  .  This  expression 
poses  considerable  difficulty  from  a  computational  standpoint.  The  prob¬ 
lems  arise  basically  from  the  facts  that  Equation  (32)  contains  a  double 

2 

integral  and  approximately  N  summations  are  required  in  its  numeri¬ 
cal  evaluation  (as  compared  to  N  summations  for  a  one -dimensional 
integral)  and,  secondly,  that  the  phase  factor  containing  (h^  + 
varies  over  a  wide  range  during  the  transverse  integration  necessitating 
the  use  of  a  large  number  of  increments  in  performing  the  summation. 
Some  effort  has  been  expended  on  obtaining  an  alternate  formulation  of 
Equation  (32)  which  is  more  amenable  to  computer  programming.  The 
results  of  this  work  appears  in  Appendix  A. 

However,  a  more  fruitful  approach  has  been  developed  utiliz¬ 
ing  both  mathematical  and  physical  modifications.  This  development  is 
presented  in  Section  IV.  A. 


23 


I 


IV.  ALTERNATE  SIGNAL  PROCESSING  TECHNIQUE 


A.  Multichannel  Signal  Processing 

As  was  stated  in  the  previous  section,  the  difficulty  in  eval¬ 
uating  the  general  sampling  function  equation  is  due  primarily  to  the 
presence  of  both  space  and  time  variables  in  the  integrands  of  the  Four¬ 
ier  transform  pairs.  In  the  analysis  that  follows,  alternate  approxima¬ 
tions  are  applied  which  allow  the  time  variable  to  be  transformed  into 
the  aperture  distribution,  and  the  transverse  variable  £  is  eliminated 
using  the  directional  characteristics  of  a  linear  array  sampling  antenna. 
The  effect  of  varying  the  filter  center  frequency  is  also  examined  by 
means  of  a  more  generalized  filter  function.  The  significance  of  these 
modifications  is  discussed  below. 


Starting  with  Equation  (8),  let 


H(w) 


where 

B  =  bandwidth 

w  =  center  frequency,  n  =  1,  2,  3,  .  .  .  ,  N. 
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For  this  condition,  Equation  (8)  becomes 


(">2  rC2|Tt>  p» 

y  *  J  d"  J  ,  ,  \  —ZTZ — HHt-t’3  *• 

y  (2ir)  Si1  c^ti)  ^-00  n 


S„(t)  = 


JL 


•I 


00  exp 

k  - 

oo 


-jk*|(vt'  -  y\)Z  +  h2  +  £2 


\  (vt'  -  r|  )2  +  h2  +  i2 


exp(jw  t')  dt' 


(33 


Since 


I 


oo  sin 


-oo 


/  co  -  co  A 

.( _ HJ 

— exp[j«(t-t«)]dw  =  7rexp[-jWn(t'-t)]  t  -  ±  <  t*  <  t  +  ± 


=  0 


for  all  other  t'  , 


Equation  (33)  becomes 


nr\2  pc2ni  p 

S(t)  =  i  \  dr,  \  d|F(4,q)  \ 

Y  Jrix  Jcl(q)  Jt 


t  +  1/2B 


1/2B 


-jk|J( 


2  2  2 
exp  l-jk^vt'  -  q)  +  h  +  £ 


N  (vt1  -  r(  )2  +  h2  +  42 


expfj(co  -u  )”]t'  exp(-jw  t)  dt'  .  (3- 

|_  o  nj  n 
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With  the  change  of  variable 

t'  =  t"  +  t 
dt'  =  dt"  , 

Equation  (34)  becomes 


can  be  approximated  by- 


in  the  exponent  and 
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Nh2  +  £2  +  (r)  -  vt)2 


(38) 


t 


in  the  denominator.  These  approximations  of  course  impose  certain 

limitations  on  the  functions  vt"  and  (vt  -  ii).  In  particular,  the  ex- 
2 

pression  (vt")  has  been  ignored  completely  in  both  the  phase  and 
amplitude  factor.  It  is  seen  from  Equation  (36)  that  this  term  has  its 
greatest  effect  at  (vt  -  rj )  =  0  and  £  =  0.  This  is  comparable  to  the 
on-axis  quadratic  phqse  error  of  a  linear  array,  where  in  this  case,  the 

y 

"linear  array"  is  the  flight  path  —  as  illustrated  in  Figure  5.  The  max- 

13 

v 

imum  dimension  of  the  "array"  is  —  ,  and  if  the  quadratic  phase  error 

B  4 

is  to  be  ignored,  this  dimension  must  be  limited. 


Figure  5.  Variations  in  R  Associated  with  the  Flight  Path 
Length  — 
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Referring  to  Figure  6, 


h2  +  (rk)  =  r2  =  <h  +  ar)2 

2  /v  f  2  2 

h  +  [zb )  =  h  +(AR)  +  2hAR  • 

2 

Where  AR  is  small  enough  to  allow  the  deletion  of  (AR)  ,  then 


2hAR 


Using  a  —  criteria  for  the  maximum  allowable  phase  error  produces 


v 


Figure  6.  Geometry  for  Determining  the  limits  on  B 
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This  forces  the  condition 


B  > 


v 

4hT 


on  B  .  For  v  =  880  feet  per  second,  h  =  1580  feet,  and  X.  =  2  feet 


B  > 


880 

*43160 


1 5.6  cycles 


which  is  compatible  with  the  conditions  of  Section  III.  A.  Since  it  has 
been  shown  that  the  range  of  (q  -  vt)  is  limited  in  the  £,  q  plane,  Equa¬ 
tions  (37)  and  (38)  are  comparable  to  the  usual  far  field  approximations. 
Under  these  conditions,  Equation  (35)  can  be  integrated  with  respect  to 
t"  to  give 


S  (t) 

y 


j  exp  (ju^t)  0^2  P 

"  - 4rrB  J  dT1  J 


c2(ri) 


d£F(£,  q)  k 


c  x  ( T1 ) 


r 


exp 


[-jktjh2  +£2  +(ti  -  vt)2 ] 


sin 


2B 


kv(  T|  -  vt) 


*Jh2  +  |2  +(r|  -  vt)2 


+  w  -  w 
o  n 


J-  .  .2 

N|h  +|  +(r)  -  vt) 


1 

2B 


kv(r|  -  vt) 


4 


h2+|2+(T1  -vt)2 


+  w  -  w 
o  n 


(39 


Equation  (39)  differs  from  Equation  (15)  in  that  the  term  (q  -  vt)  appears 

under  the  radical  and  the  term  w  -  w  is  present  in  the  — —  sampling 

on  u 

function.  The  significance  of  the  former  will  become  more  apparent  as 
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the  analysis  continues.  Regarding  the  term  it  is  apparent  that 

this  factor  shifts  the  peak  of  the  8in  U  functions  in  the  ri  direction 
ahead  or  behind  the  zero  Doppler  frequency  trace  as  is  greater  than 
or  less  than  .  Figure  7  illustrates  the  result  of  passing  the  received 
signal  through  a  multichannel  filter  bank.  By  proper  choice  of  filter 

frequencies,  the  £,  ti  plane  has  been  divided  into  strips  where  the  out- 

til  til 

put  of  the  n  filter  is  the  signal  transmitted  from  the  n  Doppler  fre¬ 
quency  trace.  The  presence  of  the  term  (r|  -  vt)  verifies  mathematically 
that  the  sampling  system  is  of  a  traveling  nature;  i.  e.  ,  with  the  excep¬ 
tion  of  the  aperture  distribution  function,  F(£,  r|),  each  factor  in  the  in¬ 
tegrand  of  Equation  (39)  is  a  traveling  function  with  velocity  v  in  the 
direction. 


Figure  7.  Sampling  Strips  Produced  by  a  Multichannel 
Filter  Bank 
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At  this  point  it  is  advantageous  to  introduce  the  effect  of  re¬ 
placing  the  omnidirectional  sampling  probe  with  a  linear  array.  Figure 
8  illustrates  the  geometry  of  a  linear  array  and  its  trace  in  the  (£»ti) 
plane.  The  array  is  oriented  at  the  point  (0,0, h)  in  the  £,ti,4  coordi¬ 
nate  system  with  its  longitudinal  axis  parallel  to  the  f  axis.  For  an 
assumed  S*n  "■  pattern,  the  argument  u  takes  the  form 


u 


irL 

X 

o 


sin  (3 


where  L  is  the  length  of  the  array  and  (3  is  the  broadside  angle  measured 
from  the  £,ti  plane. 


Figure  8.  £,  Trace  of  a  Linear  Array  Sampling  Antenna 


31 


4 


Under  these  conditions  the  weighting  effect  of  the  array  on  a 
signal  received  from  the  plane  can  be  expressed  as, 


(4( 


where  the  space  attenuation  factor 


1 


2 

+  *1 


has  been  included.  Equation  (40)  can  now  be  generalized  by  rotating  the 
trace  about  the  £  axis  through  an  angle  a  and  then  translating  the  skewed 
strip  along  the  q  axis  by  the  time  function  vt.  The  resulting  expression 


sin, 

\\ T  — 

ttL 

£  cos  a  -  (r)  -  vt)  sin  a 

X 

/\jh2  +  §2  +  (q  -  vt)2 

2‘  ,l[ 

£  cos  c  -  (ti  -  vt)  sin  a 

X 

N  h2  +  £2  +  (r|  -  vt)2 

represents  the  trace  of  a  linear  array  carried  by  the  mapping  vehicle  at 
an  angle  a  with  respect  to  the  direction  of  motion.  Under  these  condi¬ 
tions  the  expression  for  the 'filtered  signal  is 
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j  exp  (ju  t)  Pn 2  PC2(,1)  exp  ( -jkr 1 

T1 ,  '■'cJr,)  (r1) 


(41 


'1  1 


- 

>1 

r 

sin  . 

1 

kv(  T1  -  vt) 

.  sin  . 

uIj 

2B 

2  ?  0  n 

X 

\ 

h-  +  4' +(ti  -vtr 

L. 

i—N 

J 

|  cos  a  -  (q  -vt)sinc 


[_  “n!  h2  +  42  +(ti  -  vt)2 


1 

2B 


kv(ri  -  vt) 

- - - - -  +  W  -  OJ 

7TL 

4  cos  a  -  (tj  -  vt)  sin  a~| 

2  ?  2  0  n 
h‘ +4  +(h-vtr 

X 

_  x 

,2  ,  .2 
h  +4  +( t|  -  vt) 

d4dri 


where 


(r1)2  =  h2  +  42  +(ri  -  vt)2 


Figure  9  illustrates  the  limiting  effect  of  both  the  filter  function  and  the 
pattern  factor  of  the  sampling  array.  It  is  seen  that  the  array  has  re¬ 
duced  the  sampling  strip  produced  by  the  filter  to  a  sampling  increment 
with  dimensions  in  the  %  direction  comparable  to  those  in  the  r|  direc¬ 
tion.  Hence  far  field  approximations  can  now  be  made  in  the  function 


exp  (-jkr1) 
(r-)2 


with  respect  to  both  4  ar*d  ri  ■ 

In  order  to  select  a  suitable  coordinate  system  for  the  far 

field  expression,  attention  is  directed  to  the  point  (4  *  n1  +  vt)  in 

o  o 

Figure  9.  This  point  is  the  intersection  of  the  two  lines  along  which  the 
two  8in  u  terms  in  Equation  (41)  have  maximum  values.  Hence,  for 
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Figure  9.  Incremental  Sampling  Area  Isolated  by  the  Filter 
Function  and  Sampling  Array  Pattern  Factor 


given  values  of  L,  u  ,  u  ,  and  a  ,  the  coordinates  £  and  r»  ’  can  be 

on.  o  o 

computed  by  equating  each  — — —  argument  to  zero  and  then  solving  the 
two  simultaneous  equations  for  f;  and  -vt  as  follows.  From 


1 

2B 


kv(  T)  -  vt) _ 

.2  2 
+  £  +  (T)  -  vt) 


+  w 

o 


0 


and 


ttL 


£  cos  a  -  (r|  -  vt)  sin  a 
'Jh2  +  %Z  +  (ti  -  vt)2 
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(42 


(»)  -  vt)  = 


,2.  x2 

h(Wn-Wo)  .  ..2 

- __  =  (11.) 


2  2 
(kv)  -  (w  -  w  )  sec  a 
n  o 


2  2  2 
h  («  -  u  )  tan  a  - 

f  2  _  n  o _  -  t 2 

*  =  .2  .2  2  = 


(43 


(kv)  -  (to  -  w  )  sec  a 
n  o 


The  far  field  approximations  in  polar  coordinates  can  now  be  constructed 
for 

exp  ( -  jkr ') 

(r1)2 

at  the  origin  and  then  translated  to  the  point  (£o>  +  vt,  0) 

as  follows: 

4 

At  the  point  (0,  0,  0) 


e.xJ* ,  1  ~ »  exP  (  jkR)  exp  sjn  Q  cos  (jj  +  r|  sin  <jT|1 
(r'j  R  1  JJ 


where  R  is  a  constant  slant  range  and  0  and  <)>  are  the  familiar  azimuth 
and  coelevation  coordinates  illustrated  in  Figure  10. 
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Figure  10.  Polar  Coordinates  for  the  Far  Field  Approximations 
Referred  to  the  Origin 


At  the  point  (g^,  +  vt,  0) 


exp  (  -jkr  )  ^  exp  (  -jkR)  exp  fjk  sin  ef(  |  -  £  )  cos  <j>  +  Qi  -(n^  +  vt)J  sin  <p"l  ^ 
(r')  R  V 


where  now 
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2  .  2  ,2 
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The  result  of  this  translation  is  illustrated  in  Figure  11.  Equation  (41) 
now  becomes 


(44) 
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Figure  11.  Polar  Coordinates  for  the  Far  Field  Approximations 
Referred  to  the  Point  (£q,  T|^+vt,  0) 
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Although  mathematically  unwieldy,  Equation  (45)  yields  considerable  in¬ 
formation  as  tc  the  nature  of  the  sampling  system  output.  From  previous 
experience  with  the  analysis  of  illuminated  apertures.  Equation  (45)  is 
recognized  as  the  far  field  diffraction  integral  of  an  irregular  shaped 
aperture  with  a  tapered  distribution.  The  weighting  functions,  and 

,  clearly  have  a  tapering  effect  on  the  distribution  function  F(4»  "n )  , 
and  for  a  chosen  attenuation  criteria  establish  the  boundaries  of  the  iso¬ 
lated  area  in  the  4>  *1  plane.  It  is  again  emphasized  that  the  piesence 
of  the  term  (q  -  vt)  establishes  the  traveling  nature  of  the  integrand; 

i.e.  ,  the  isolated  area  centered  at  the  point  (4  .  V  +  vt,  0)  ,  the  slant 

o  o 

range  R  ,  and  the  vehicle  position  vt  are  all  fixed  with  respect  to  each 
other  and  travel  at  a  constant  velocity  v  in  the  q  direction. 

Because  of  the  time  dependence  of  the  sampling  function,  the 
limits  of  the  integral  and  hence  the  distribution  function,  F(4>q)  ,  are 
also  functions  of  time.  This  can  be  illustrated  by  making  the  approxi¬ 
mations  in  the  arguments  of  the  Sm  U  functions,. 


„Z  ,  .2 

+  4  +  (q  -  vt) 


R 


(4 


which  is  valid  for  small  variations  in  4  and  q  about 


4  and  q  1  . 
o  o 


Now 
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Figure  12.  Incremental  Sampling  Area  Using  a  Constant 
Value  of  R  in  the  Filter  and  Array  Functions 
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Using  the  3  db  attenuation  criteria,  the  limits  on  £  and  r| 
are  computed  from 


1  kv(  T|  -  vt) 

2B  L  R 


+  (jj 

o 


±1.39 


and 


ttL  %  cos  a  -  (ri  -  vt)  sin  a 

\  L  R 


±1.39 


From  which  follows 


vt  =  ±  —  2.78B  +  (lo 
kv  kv  n 
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kv 
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and 
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±  1.39  — p  sec  a  +  ( t)  —  vt)  tan  a 
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Under  these  conditions,  the  limits  of  integration  in  Equation  (42)  are 


where 


tl .  =  vt  +  -n  *  -  A  2.78  B 

1  o  kv 

=  vt  +  T|  *  +  -£-  2.78  B 

2  o  kv 

It  is  seen  that  the  limits  on  £  are  a  linear  function  of  t)  and  t,  and  r| 
varies  symmetrically  about  a  particular  value,  (vt  +  r| 1  )  ,  which  is  a 
function  of  the  time  variable.  Hence  the  distribution  function  F(£,t|)  is 
necessarily  a  function  of  t  because  of  the  nature  of  the  limits  of  integra¬ 
tion.  This  result  can  be  obtained  in  a  more  straightforward  manner  by 
the  change  of  variable 


T|  =  T)  1  +  Vt  .  (5 

Substituting  the  conditions  of  Equations  (48)  and  (55)  into  Equation  (45) 
produces 


j  exp[j(oj  t  -  kR jl  0^2  PC2^r|  ^ 

g(t)  =  - ■■  -  \  \  F(  1!’  +vt) 

4-rrBR‘i  ■  Uc'(t,') 

•  exp^jk  sin  6  jj  £  -  4q)  cos  <j>  +  (r|  '  -  r|^)  sin  W ^  d£dr|  1 
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It  is  seen  that  the  change  of  variable  of  Equation  (55)  has  essentially 
transformed  the  time  variable,  t,  into  the  aperture  distribution  func¬ 
tion,  while  at  the  same  time  removing  it  from  the  remaining  factors 
and  limits  of  the  diffraction  integral.  This  is  in  complete  agreement 
with  the  foregoing  discussion;  i.e.  ,  at  any  instant  of  time  the  output  of 
the  filter  represents  the  weighted  aperture  distribution,  F(|,  q1  +  vt)  , 
integrated  over  the  limits. 
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Figure  13.  Geometry  of  the  Moving  Reference  Frame,  (4,  q',  l,), 
and  the  Fixed  Frame  (4>  r| ,  £) 


Equations  ( 56)  through  (  59)  can  be  interpreted  as  the  far  field  expression 

for  a  radiating  aperture  referred  to  as  a  moving  reference  frame,  i.  e., 

the  coordinate  frame  of  the  sampling  system.  In  the  4>  q '  plane,  the 

incremental  area  is  a  parallelogram  centered  at  the  point  4  =  and 

q  '  =  q^  as  illustrated  in  Figure  13.  Equations  (42)  and  (43)  give  the 

values  of  q  '  and  4  in  terms  of  the  parameters  h,  w  ,  w  ,  v,  and  a  . 
o  o  on 

Hence,  it  follows  that  for  fixed  values  of  h,  w  ,  v,  and  a,  variations 

o 

in  the  parameter  w  produce  a  shift  in  the  position  of  the  point  (4^,q^) 

along  the  line  4  =  tj  1  tan  a  .  For  the  case  of  the  multichannel  filterbank 

discussed  previously,  the  trace  of  the  linear  array  is  divided  up  into  N 

th 

sampling  increments.  The  output  of  the  n  filter  is  in  the  form  of 
Equation  (53)  where  the  aperture  distribution  is  now  clearly  a  function 
of  time.  The  output  of  the  n^  filter  can  now  be  expressed  as 
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itL 
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exp  ^jk  sin  0  j(  £  -  cos  <f>  +  (V  “il^)  sin  <t>jj  F(l>  *n  '  +vt)  W^W2  (6 


where  W  and  W  differ  from  Equations  (57)  and  (58)  only  in  that  R  is 

1  Z  th 

replaced  by  R  to  denote  the  n  filter  output. 

Figure  14  illustrates  the  geometry  of  a  multichannel  mapping 

system.  It  is  apparent  that  at  any  instant  of  time,  the  field  distribution 
th 

over  the  n  incremental  sampling  area  corresponds  to  that  portion  of 

the  q  plane  which  lies  within  the  boundaries  of  the  sampling  increment. 

Hence  it  would  appear  that  each  filter  output  provides  a  continuous  strip 

mapping  of  the  radiating  aperture  as  the  mapping  vehicle  traverses  the 

q  plane.  Continuing  along  this  line  of  reasoning,  it  would  follow  that 

after  suitable  correction  for  the  weighting  factors,  variation  in  R  ,  etc., 

n 

the  combined  output  of  the  multichannel  filter  bank  would  result  in  a 
complete  mapping  of  the  aperture  distribution  in  the  §,  t|  plane.  Since 
the  aperture  distribution  is  continuously  averaged  over  the  dimensions 
of  the  parallelogram -shaped  sampling  increment,  the  accuracy  of  the 
measured  distribution  would  of  course  depend  upon  the  attainable  reso¬ 
lution.  Unfortunately  other  factors,  which  at  this  point  are  not  apparent, 
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Figure  14.  Sampling  Areas  of  a  Multichannel  Mapping  System 

enter  into  the  approach  outlined  above  which  greatly  restrict  the  posi¬ 
tion  of  the  sampling  increment. 

Equation  (60)  can  be  analyzed  further  on  a  qualitative  basis 
by  making  the  following  simplifications:  If  the  aperture  distribution 
F(£,  r\)  is  assumed  to  be  essentially  constant  over  the  limits  of  integra¬ 
tion,  the  term  F(£,  q  1  +  vt)  can  be  approximated  at  each  instant 

of  time  by  an  average  value  A(t) .  This  amounts  to  assuming  a  uniform 
distribution  over  the  sampling  increments  which  varies  uniformly  in 
intensity  with  t.  Under  this  substitution  it  is  possible  to  perform  the 
integration  of  Equation  (60);  however,  it  is  recognized  that  the  tapering 
effect  of  the  weighting  functions  and  would  tend  to  broaden  the 
beam  and  decrease  the  sidelobe  level  of  the  resulting  diffraction  pattern 
Under  the  assumptions  stated,  Equation  (60)  becomes 


45 


gn(t)  =  C  A(t) 


R 

+  (2.78B) 


R 

rr(2.78B) 


sin  ©  ( r| 


ti  1 )  sin 
n 


■I 


R  X. 

■n  1  tan  a  +  1.39  — 7-  sec  a 
itL 


■p  '  tan  a 


R  \ 

1.39  — ^7-  sec  a 

IT  Li 


exp^jk  sin  0(£  -  %  )  cos  ((Tj  d§ 


(61 


where 


j  expfj 

ICO  t  -  kR  n 

L 

O  n  J 

4BirR2 

n 


From  Figure  11  the  following  relations  are  apparent 
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Under  these  conditions  Equation  (61)  becomes 
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Performing  the  integration  with  respect  to  £  produces 
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Integrating  with  respect  to  T) 1  produces 
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which  can  be  written  as 
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The  effect  of  and  £  on  the  received  signal  is  immediately  apparent 

from  Equation  (62).  It  is  clear  that  the  signal  g  (t)  is  maximum  for 

n 

r|^  =  £n  =  0  independent  of  angle  a  .  It  is  also  clear  that  as  the  radial 
distance 


P  = 


sin  u 

increases,  the  received  signal  is  attenuated  by  the  two  -  terms 

u 

through  a  series  of  nulls  and  relative  maximums.  To  illustrate  this 
effect,  Equation  (62)  has  been  computed  in  the  form  of 


«p>  -  (^1  • 

This  corresponds  to  the  following  parameter  values 


(6 
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_1_ 

L 


where 


B  =  16  cycles 
v  =  1000  feet/second 
L  =  41.6  feet 
y  =  0.056  . 

Figure  15  is  a  curve  of  Equation  (63)  where  the  radius  dis- 

o 

tance  p  is  now  measured  along  the  37  diagonal  passing  through  the 

point  (i  ,  q  ' )  and  the  point  (0,  0)  in  the  £,  q  '  plane.  It  can  be  concluded 
n  n 

from  these  results  that  the  sampling  increment  must  be  limited  to  a 
position  directly  under  the  probe.  The  signal  from  any  other  increment 
centered  outside  this  region  will  be  attenuated  in  a  manner  illustrated 
by  Figure  15  which  would  clearly  make  the  accuracy  of  the  processed 
signal  questionable  and  in  some  cases  nonexistent.  Under  these  restric¬ 
tions,  Equation  (60)  can  be  greatly  simplified  by  performing  the  follow¬ 
ing  substitutions.  The  center  frequency  of  the  sampling  filter  is  taken 

to  be  the  transmit  frequency  w  .  No  apparent  advantage  is  gained  by 

o 

turning  the  sampling  array,  hence  the  skew  angle  a  can  be  equal  to  zero. 

Finally,  the  coordinates  6  and  ri '  are  now  £  =  £  =0  and  n '  =  ri '  = 

n  n  no  no 

0  .  Under  these  substitutions,  Equation  (60)  becomes 


gQ(t)  = 


j  exp  [j(«  t  -  khj] 


BTTh 


Si  [l.39]  Si  [l.39]  A(t) 


(6< 


49 


Figure  15.  Signal  Strength  Versus  Sampling  Increment  Displacement 


where  Si  is  the  tabulated  sine  integral.  Equation  (64)  reduces  to 

A(t)  =  j2Bh2  expp(uQt  -  kh)J  gQ(t) 

where  A(t)  is  the  value  of  the  aperture  distribution,  F(£,r|),  averaged 
over  a  rectangular  sampling  increment  centered  at  the  point  4=0, 
r]  =  vt .  Figure  16  illustrates  the  geometry  of  this  final  configuration. 


(65) 
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Figure  16.  Principal  Plane  Mapping  by  Means  of  a  Single 
Incremental  Sampling  Area 


While  the  feasibility  of  mapping  the  complete  aperture  in  one  flight 
transit  has  proven  to  be  unrealistic,  it  is  apparent  that  conventional 
plane  cuts  can  be  obtained  by  making  several  passes  over  the  radiat¬ 
ing  aperture.  The  signal  obtained  from  each  transit  would  be  a  con¬ 
tinuous  mapping  of  the  strip  traced  out  by  the  sampling  increment.  As 
was  stated  before,  the  accuracy  of  the  measured  distribution  depends 
mainly  upon  the  size  of  the  rectangular  sampling  increment;  i.e.,  the 
finer  the  resolution  the  more  nearly  the  functioA  A(t)  represents  the 
distribution  function  at  the  point  (0,vt)  .  Also,  the  limitations  on  the 
position  of  the  sampling  increment  have  been  shown  to  be  a  direct  result 
of  the  size  of  sampling  increment. 
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This  last  observation  can  be  further  illustrated  by  a  simple 
application  of  diffraction  theory  to  the  conditions  of  the  isolated  samp¬ 
ling  increment.  Consider  a  1000-foot  circular  aperture  which  is  to  be 
measured  by  a  vehicle  traveling  at  an  altitude  of  1000  feet  over  the  cen¬ 
ter  of  the  radiating  plane.  Each  isolated  portion  cf  the  aperture  can  be 
thought  of  as  exhibiting  a  diffraction  pattern  which  is  a  function  of  the 
increment  size.  If  the  mapping  vehicle  is  contained  within  the  half  power 
beamwidth  of  the  most  distant  increment,  it  is  clear  that  the  beamwidth 
of  the  increment  must  be  of  the  order  of  2  tan  *  ~~  =  53°.  This  cor¬ 
responds  to  a  radiating  aperture  with  maximum  dimensions  of  less  than 
a  wavelength. 

Several  factors  affect  the  resolving  power  cf  the  sampling 
system.  Increasing  the  altitude  of  the  probe  decreases  the  off-axis  angle, 
but  this  effect  is  nullified  by  the  increased  area  subtended  by  the  samp¬ 
ling  antenna  beamwidth.  Increasing  the  velocity  of  the  sampling  vehicle 
would  produce  a  greater  Doppler  shift  and  hence  a  greater  filter  resolu¬ 
tion  for  a  given  bandwidth,  but  this  occurs  in  one  dimension  only,  and  it 

v 

has  been  shown  that  the  ratio  —  is  subject  tc  a  maximum  critical  value. 

B  J 

One  possible  approach  to  this  problem  is  synthetic  aperture 
type  signal  processing.  This  is  analyzed  in  some  detail  in  the  following 
section. 
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B.  Synthetic  Aperture  Processing 


One  possible  approach  to  obtaining  greater  sampling  resolu¬ 
tion  is  the  application  of  synthetic  aperture  type  processing  to  the  signal 
received  from  the  sampling  area.  In  recent  years  this  system  has  been 
developed  in  connection  with  ground-mapping  radar  to  such  a  degree  that 
azimuth  resolutions  (independent  of  range  and  wavelength)  of  the  order 
of  a  few  feet  are  obtainable.  Since  most  of  the  material  available  on  this 
subject  is  still  classified,  only  a  brief  discussion  will  be  presented  here. 

In  ground-mapping  radar  applications,  the  transmitting 
antenna  used  in  a  typical  synthetic  aperture  system  is  "side  looking", 
relatively  narrow  beam,  and  of  physical  dimensions  compatible  with 
airborne  requirements.  Range  resolution  is  obtained  by  conventional 
pulsing  techniques  while  azimuth  resolution  is  obtained  by  recording  a 
series  of  returns  from  each  point  target  and  then  processing  the  recorded 
data  such  that  the  physical  antenna  becomes  an  element  of  a  large  syn¬ 
thetic  array.  The  number  of  elements  in  the  array  is  equal  to  the  num¬ 
ber  of  pulses  recorded  from  each  point  target.  Figure  17  illustrates 
the  geometry  of  a  typical  airborne  ground-mapping  radar. 

It  is  apparent  that  the.  number  of  return  pulses  and  hence  the 
length  of  the  synthetic  array  is  a  function  of  the  pulse  repetition  rate  and 
the  length  of  time  that  the  point  target  at  P  lies  within  the  specular  area 
of  the  antenna  beam.  Hence,  flight  velocity  and  beamwidth  are  also 
factors  governing  the  length  of  the  synthetic  array.  Broadly  speaking, 
the  length  of  the  synthetic  array  will  be  comparable  to  the  longitudinal 
dimension  of  the  specular  trace  and  will  have  the  conventional  resolving 
power  of  a  physical  antenna  of  that  length. 
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►  V 


Figure  17.  Geometry  of  a  Side -Looking  Synthetic  Aperture 
Processing  Ground -Mapping  Radar 


In  developing  a  suitable  processing  system  it  can  be  assumed 
that  the  coherently  pulsed  signal  is  continuous.  Under  this  condition, 
the  returned  signal  from  a  particular  point  target  is  continuous  with  its 
frequency  modified  by  a  time  and  range  dependent  Doppler  effect.  There¬ 
fore  each  point  target  generates  its  own  continuous  and  unique  Doppler 
history  as  it  passes  through  the  specular  trace  of  the  physical  antenna. 
For  example,  in  the  case  of  a  broadside  looking  antenna  (i.e.  ,  main 
beam  normal  to  the  direction  of  flight),  the  frequency  of  the  received 
signal  changes  from  a  higher  to  a  lower  value  than  that  of  the  transmitted 
frequency  as  the  target  enters  and  leaves  the  specular  trace.  Further, 
when  the  target  is  exactly  at  the  broadside  position,  the  transmit  and 
receive  frequencies  are  the  same.  This  configuration  has  been  used  in 
connection  with  a  coherent  reference  oscillator  and  a  pulse  integration 
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system  in  which  the  complete  return  from  each  target  is  compressed 

about  the  "zero  beat"  point  produced  by  beating  the  time -varying  return 

6 

signal  against  the  stable  local  oscillator.  The  zero  beat  point  and 
hence  the  center  of  the  compressed  return  signal  serves  to  identify  the 
azimuthal  position  of  the  point  target  with  respect  to  the  mapping  vehicle . 

Processing  techniques  have  been  developed  using  vector 

addition,  cross  correlation,  and  filtering  approaches.  However,  in  a 

recent  publication  it  has  been  shown  that  these  three  approaches  are 
7 

equivalent,  hence  in  analyzing  a  particular  aspect  of  a  given  system, 
the  most  applicable  approach  should  be  used. 

Although  grossly  over-simplified,  the  foregoing  discussion 
contains  sufficient  information  for  deducing  certain  conclusions  regard¬ 
ing  the  mapping  of  radiating  apertures  via  synthetic  aperture  processing 
techniques.  As  in  any  pulsed  radar,  the  operation  of  a  synthetic  aper¬ 
ture  system  consists  of  periods  of  transmission,  during  which  time  wave 
trains  of  finite  length  are  launched  towards  the  target  area,  followed  by 
periods  of  reception  during  which  time  the  reflected  portions  of  the  trans¬ 
mitted  wave  are  received.  Hence  during  the  period  of  reception,  the 
target  area  is  analogous  to  a  two-dimensional  array  whose  elements 
consist  of  all  the  point  targets  which  lie  within  the  specular  trace  of  the 
physical  antenna.  Therefore,  it  would  appear  that  the  transition  from 
ground  mapping  to  aperture  mapping  is  a  logical  and  straightforward 
step.  However,  further  investigation  into  the  nature  of  the  quasi-point 
sources  on  which  the  analyses  of  the  two  "radiating  apertures"  are  based, 
brings  out  complications  which  appear  to  make  the  two  systems  incom¬ 
patible. 
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Continuing  the  analysis  of  the  specular  trace  on  the  basis  of 
a  two-dimensional  array,  it  is  apparent  that  the  energy  .distribution  on 
each  element  (i.e.  ,  point  target)  is  dependent  upon  the  reflection  coef¬ 
ficient  and  the  position  of  the  point  target  with  respect  to  the  illuminating 
and  receiving  antenna.  Hence  the  target  area  as  an  array  is  made  up  of 
an  arbitrary  number  of  elements  with  distributions  essentially  random 
in  both  magnitude  and  phase.  Clearly  this  is  a  necessary  condition  for 
the  operation  of  synthetic  aperture  type  processing  (i.e.  ,  each  target 
must  have  a  unique  Doppler  history  during  its  passage  through  the  spec¬ 
ular  trace).  For  example,  consider  a  target  area  made  up  of  a  large 
number  of  Identical,  closely-spaced  point  targets.  Under  these  condi¬ 
tions  the  reflected  energy  is  no  longer  made  up  of  discrete  incoherent 
wavelets.  Instead,  as  the  separation  between  targets  is  decreased,  the 
reflected  wave  trains  from  adjacent  targets  begin  to  overlap  in  a  well- 
defined  manner.  The  result  is  threefold. 

1.  The  overlapping  of  the  returned  wave  trains  makes  pulse 
ranging  impossible. 

2.  The  Doppler  frequency  histories  cf  adjacent  targets  begin 
to  "look"  alike. 

3.  Reinforcement  and  cancellation  resulting  from  the  coher¬ 
ence  of  the  reflected  wavelets  produce  the  usual  lobed 
pattern  associated  with  any  aperture  which  is  illuminated 
in  an  orderly  manner. 

Any  of  these  three  effects  would  make  the  synthetic  processing  system 
inapplicable.  However,  for  an  actual  illuminated  aperture,  the  third 
effect  as  discussed  in  the  preceding  section  is  probably  the  most  deci¬ 
sive.  Hence,  it  appears  that  the  coherent  nature  of  the  quasi-point  dis¬ 
tribution  of  a  radiating  aperture  is  incompatible  with  the  requirements 
of  synthetic  aperture  processing. 


56 


Having  established  the  limitations  of  the  foregoing  analysis 
as  regards  the  number  and  position  of  sampling  areas,  any  further  dis¬ 
cussion  will  be  restricted  to  the  conditions  of  the  area  directly  under  the 
probe.  In  general,  the  distribution  over  an  isolated  portion  of  an  arbi¬ 
trarily  illuminated  aperture  will  be  nonlinear  and  asymmetrical  in  both 
amplitude  and  phase.  The  added  weighting  effect  cf  the  filter  function 
and  sampling  array  pattern  produces  an  aperture  distribution  that  is 
extremely  complex.  Fortunately,  only  a  qualitative  analysis  is  needed 

and  this  can  be  accomplished  with  a  sufficient  degree  of  accuracy  from 

8 

line  source  considerations.  As  reported  by  R.  C.  Spencer,  the  diffrac¬ 
tion  field  of  a  line  source  can  be  analyzed  in  terms  of  the  line  distribu¬ 
tion  moments.  For  an  arbitrary  distribution,  f  (x',  the  normalized 
diffraction  integral  is, 
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Equation  (66)  can  be  expanded  into  its  infinite  series  with  the  result 
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is  referred  to  as  the  k  moment  of  f(x)  .  The  power  pattern  P(u) 
follows  from  the  sum  of  the  squares  of  the  real  and  imaginary  series  of 
Equation  (67), 
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Multiplicative  constants  have  been  deleted  in  Equations  (66)  through  (68). 


Several  important  factors  are  immediately  apparent  from 
Equations  (67)  and  (68).  First,  it  follows  from  the  definition  of  p.^  that 
even  and  odd  distribution  functions  produce  respectively,  pure  real  and 
pure  imaginary  diffraction  field  series,  which  are  again  respectively 
even  and  odd  functions  of  the  variable  u.  In  the  general  case,  both 
real  and  imaginary  components  are  present,  however.  Equation  (68) 
indicates  that  the  power  pattern  is  still  symmetric  about  u  =  0  and  is  a 
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principal  maximum  there.  Thus,  the  only  asymmetry  resulting  from  an 

arbitrary  amplitude  distribution  is  the  sign  reversal  in  the  imaginary 

component  of  the  field  intensity  which  accompanies  the  sign  reversed  of 

the  u  variable.  The  effect  of  the  linear  phase  distribution  can  be  seen 
4 

from  the  relation 


it  a.  ,  _ 

u  =  —  sin  0  -  (3 
\ 

For  u  =  0  ,  ~  sin  9  =  (3 ,  hence  the  principal  maximum  or  direction  of 
the  main  lobe  is  skewed  by  an  angle  (3  . 

From  these  results,  the  following  assumptions  can  be  made 
regarding  large  apertures  with  arbitrary  amplitude  distributions  and 
moderate  but  continuous  phase  variations.  A  small  isolated  area  cen¬ 
tered  about  a  point  in  the  plane  of  the  aperture  exhibits  a  symmetrical 
power  pattern  with  the  principal  lobe  axis  normal  to  the  phase  front  at 
the  point.  The  on-axis  field  intensity  is  proportional  to  the  zero  distri¬ 
bution  moment  which  is  simply  the  aperture  distribution  integrated  over 
the  incremental  sampling  area,  and  the  extraction  of  relative  phase  and 
amplitude  from  a  signal  of  this  type  can  be  performed  in  a  conventional 
manner. 
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V. 


ERROR  ANALYSIS 


The  errors  associated  with  performing  pattern  measurements  with 
the  Transverse  Doppler  Pattern  Measurement  Technique  fall  logically 
into  two  categories:  one  includes  the  inaccuracies  which  are  inherent  in 
the  technique:  the  other  consists  of  those  errors  which  arise  from  devia¬ 
tions  in  the  assumed  values  of  system  parameters. 

The  former  are,  in  general,  functions  of  the  t  /stem  aperture  re¬ 
solution,  where  the  resolution  is  dependent  on  the  composite  parameter 
kv£i 

longitudinal  dimension,  and  on  the  directional  characteristics 
of  the  sampling  probe  in  the  transverse  dimension.  Since  the  resolution 
increment  is  of  finite  dimension  and  the  measuring  technique  involves 
essentially  an  averaging  of  the  response  received  from  the  increment, 
the  ability  of  the  system  to  resolve  local  variations  from  the  average  is 
limited.  In  the  case  of  the  line  source,  the  response  is  averaged  only 
over  the  longitudinal  dimension  of  the  sampling  increment.  For  a  two- 
dimensional  aperture,  sampled  by  an  omnidirectional  probe,  the  response 
is  averaged  over  the  transverse  dimension  of  the  aperture  as  well,  and 
the  filtered  signal  is  weighted  by  both  the  aperture  distribution  in  the 
transverse  direction  and  the  transverse  extent  of  the  aperture. 

However,  if  a  linear  array  is  used  for  the  sampling  antenna,  the 
transverse  dimension  over  which  the  signal  is  averaged  is  reduced  to  a 
length  proportional  to  the  beamwidth  of  the  linear  array.  Under  this 
condition,  the  weighting  effect  of  the  aperture  distribution  in  the  trans¬ 
verse  dimension  is  lessened;  however,  the  effect  of  the  pattern  factor  of 
the  linear  array  must  be  accounted  for  in  the  final  result. 
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To  determine  the  effect  of  errors  attributable  to  system  limitations 
in  the  measurement  of  the  aperture  distribution,  far  field  patterns  were 
computed  from  the  "measured"  aperture  distributions  of  a  line  source 
500X.  in  length  for  several  values  of  system  parameters  and  aperture  dis¬ 
tributions.  These  cases  are  plotted  and  discussed  on  pages  41  through 
57  of  Scientific  Report  No.  1."*  In  general,  the  patterns  computed  from 
the  measured  distributions  are  in  good  agreement  with  the  exact  patterns 
in  the  main  beam  region,  but  are  optimistic  compared  to  the  exact  pat¬ 
tern  outside  of  this  region. 

Errors  which  are  the  result  of  variations  in  system  parameters, 
e.g.,  height,  velocity,  operating  frequency,  etc.,  during  the  measure¬ 
ment  interval  are  somewhat  more  difficult  to  analyze  quantitatively. 
However,  what  qualitative  effects  these  deviations  will  have  are  readily 
deduced  from  a  consideration  of  Equation  ( 15).  Basically  the  system 
parameters  which  may  be  expected  to  vary  appear  in  three  places  in 
this  expression:  in  the  amplitude  factor,  in  the  phase  factor,  and  in  the 
resolution  of  the  sampling  function. 

Since  the  time  required  for  taking  the  measurement  data  (or  aper¬ 
ture  fly-by  time)  will  be  of  the  order  of  two  seconds,  variations  in  the 

transmitted  frequency  during  this  time  will  be  very  small,  perhaps  of 

8  9 

the  order  of  one  part  in  10  at  UHF.  The  frequency  dependence  appears 

as  the  parameter  k  in  all  three  factors  of  the  equation,  and  since 

Ak  =  k  ^  ,  variations  of  thi3  order  in  frequency  have  little  effect  on 

system  accuracy.  Even  in  the  phase  factor  a  variation  in  k  of  approxi- 

3 

mately  two  parts  in  10  are  required  for  it  radians  of  phase  change  with 
h  =  1500  feet.  An  additional  condition  on  frequency  stability  has  been 
assumed  in  the  filter  function  of  Equation  (10),  i.  e. 
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H2(o>) 


(] 


However,  since  the  band  of  frequencies  about  ,  is  essentially  contin¬ 
uous  (except  near  the  ends  of  the  aperture)  because  of  the  Doppler  con¬ 
tributions,  a  moderate  change  in  corresponds  to  sampling  a  strip  of 
the  aperture  which  is  located  either  slightly  ahead  or  slightly  behind  the 
probe  vehicle  rather  than  directly  beneath  it.  For  example,  at  a  velocity 

g 

of  880  ft/ sec,  a  variation  in  frequency  of  one  part  in  10  causes  an  angu¬ 
lar  shift  of  the  sampling  region  of  approximately  one  percent. 

Probe  height,  h  ,  appears  in  both  the  phase  factor  and  the  argu¬ 
ment  of  the  sampling  function.  In  the  sampling  function,  the  effect  of 

changes  in  h  appear  as  variations  in  the  resolution  of  the  aperture.  As 

5 

can  be  seen  from  Figures  5b,  20,  and  21  of  Scientific  Report  No.  1, 
changes  in  the  composite  parameter,  K  =  kva/2Bh,  over  a  large  range 
(2  to  1)  have  a  relatively  small  effect  on  the  accuracy  of  the  aperture 
distribution  measurement  and  the  corresponding  far  field  pattern  once  a 
certain  minimum  value  of  the  parameter  K  is  exceeded.  In  view  of  this 
fact,  it  appears  that  small  changes  in  the  values  of  any  of  the  individual 
factors  of  K  will  not  limit  the  accuracy  of  the  measurement  technique. 

The  height,  h,  appears  also  in  the  phase  factor.  The  phase  may 
be  written  as 


k*Jh2  +  |2 


kh^ 
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and  varies  essentially  as  kh .  Since  h  is  of  the  order  of  one  to  two 
thousand  feet,  very  small  changes  in  h  percentage -wise  can  cause  the 
phase  to  shift  appreciably.  For  a  linear  change  in  h,  the  effect  is  that 
of  a  linear  phase  shift  across  the  aperture  and  can  cause  difficulty. 
Therefore,  it  appears  that  it  will  be  necessary  to  limit  the  variation  in 
h  to  the  order  of  a  wavelength  during  the  measurement  period.  How¬ 
ever,  in  the  event  that  large  flight  path  deviations,  in  terms  of  wave¬ 
length,  are  unavoidable,  it  is  possible  to  detect  such  deviations  by  means 
of  a  Doppler  and/or  inertial  stabilization  system  such  as  those  used  in 
conjunction  with  ground-mapping  radars.  Such  a  system  provides  con¬ 
tinuous  monitoring  of  pitch,  roll,  altitude,  and  velocity;  hence  the  out¬ 
put  signal  can  be  used  to  compensate  for  any  deviations  in  flight  path  by 
electronically  scanning  the  sampling  antenna,  adjusting  the  filter  fre¬ 
quencies,  and  by  varying  the  phase  of  the  received  signal  in  accordance 
with  incremental  changes  in  altitude. 

As  was  pointed  out  in  Section  III,  the  presence  of  a  large  linear 
phase  distribution  in  the  radiating  aperture  necessitates  a  modification 
of  the  mapping  technique.  Fortunately  the  solution  to  this  problem  is 
straightforward  and  is  contained  in  the  general  expression  for  the  Doppler 
frequency  shift. 


Referring  to  the  geometry  of  Figure  lb,  the  frequency  of  the  signal 
received  from  a  point  in  the  £,r\  plane  by  the  probe  at  P  is  given  by 
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where  c  is  the  velocity  of  light  and  the  term 


u 


o  dr 
c  dt 


is  the  Doppler  frequency  shift. 

It  has  been  shown  that  the  radial  distance,  r  ,  connecting  P  and 
the  point  (£,  r|)  is  of  the  form 


r  = 


=  ^(vt  -  r\)Z  +  h2  +  £Z 


(7 


Therefore,  the  Doppler  frequency  equation  is, 
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It  should  be  mentioned  that  the  proper  sign  has  been  chosen  for  the  radi¬ 
cal  of  Equation  (72)  so  that  the  Doppler  shift  is  positive  for  ti  >  0  and 
negative  for  r)  <  0  . 

For  constant  Doppler  frequencies  =  w.  ,  i  =  1,  2,  3,  .  .  .  , 
Equation  (72)  becomes 
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which  reduces  to 
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as  is  the  case  here.  Equation  (74)  represents  a  family  of  hyperbolas  with 
foci  shifted  along  the  t|  axis  by  the  linear  function  of  time,  vt. 

As  was  stated  earlier,  the  physical  significance  of  this  is  that  there 
is  associated  with  a  bank  of  filters  at  the  probe  P  ,  a  grid  of  hyperbolic 
contours  in  the  §,  t)  plane  which  travels  in  the  same  direction  and  at  the 
same  velocity  as  the  probe. 

The  effect  of  a  large  linear  phase  deviation  will  now  be  considered. 
The  geometry  of  Figure  18  corresponds  to  that  of  Figure  lb  with  the  ex¬ 
ception  of  the  tilted  plane  W  which  represents  a  wavefront  emanating 
from  an  aperture  in  the  £»  q  plane  which  has  been  electronically  scanned 
through  an  angle  0  in  the  £,,q  plane.  From  the  geometry,  the  slant 

I 

range  r  connecting  the  probe  at  P  with  a  general  point  in  W  is 

r  -  \  +  (r\  -  vt)2  +  (h  -l,)2 
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Figure  18.  Plane  Wavefront  of  an  Electronically -Scanned  Aperture 


where 


Therefore 


—  =  tan  0 
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dr 
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£2  +  (t)  -  vt)2  +  (h  -  r)  tan  8)2 


Hence  the  frequency  of  the  signal  at  P  is 
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2 

Equation  (75)  differs  from  Equation  (72)  in  that  the  term  (h  -  q  tan  0) 

2 

appears  under  the  radical  instead  of  h  .  Because  of  this,  Equation  (72) 
does  not  represent  a  traveling  grid  of  fixed  contours  such  as  was  devel¬ 
oped  from  Equation  (72).  By  following  the  procedure  used  to  obtain 
Equation  (74),  it  can  be  shown  that  the  resulting  traveling  grid  is  made  up 
of  contours  in  the  £,  q  plane  which  change  shape  as  a  function  of  time. 
However,  one  important  exception  exists  which  is  more  readily  seen  if 
the  T|  plane  is  rotated  into  the  plane  of  the  phase  front.  Figure  19 
illustrates  the  results  of  this  rotation.  The  expression  for  r  is  now 


I  2  2  2 

r  =  +  (q  -  vt  cos  0)  +  (h  -  vt  sin  0) 
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Figure  19.  Plane  Wavefront  of  an  Electronically-Scanned 
Aperture  Rotated  into  the  £,t|  Plane 
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2  2  2 
+  ( r|  -  vt  cos  0)  +  (h  -  vt  sin  0) 


It  can  be  seen  from  Equation  {76)  that  the  frequency  of  the  signal  received 
from  the  point  directly  under  P  is  independent  of  t  for  if  the  §,  coor¬ 
dinates  of  P  (i.  e.  ,  |  =  0,  11=  vt  cos  0)  are  substituted  into  the  equation, 
the  result  is 
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Hence  it  can  be  concluded  that  if  the  slant  range  r  is  at  all  times  normal 
to  the  wavefront,  the  frequency  of  the  signal  received  from  the  incre¬ 
mental  area  will  be  constant  and  will  differ  from  the  fundamental  frequency 

co  by  an  amount  proportional  to  the  component  of  velocity  in  the  direction 
o 

of  r  .  Referring  to  Figure  20  it  can  be  seen  that  this  amounts  to  shifting 
the  slant  range  r  ahead  of  the  vertical  by  an  amount  equal  to  the  scan 
angle  0  . 

An  alternate  approach  to  this  analysis  is  the  consideration  of  the 
time  rate  of  change  of  phase  produced  by  the  motion  of  the  probe  with  re¬ 
spect  to  a  constant  phase  front.  Referring  to  Figure  20,  if  0  is  fixed  and 


68 


r  is  extended  to  point  P'  on  the  q  axis,  it  is  clear  that  the  length  PP1 
does  not  change  as  P  traverses  a  path  parallel  to  the  q  axis.  Hence 
the  Doppler  shift  is  zero,  however,*  the  change  in  phase  which  occurs  at 
the  moving  point  P1  is  not  zero  and  can  be  expressed  as 


d<fr 

dt 


d  2ir 

=  _  r  tl  6 

L  o  J 


r|  tan  sj  ii 


where  |r —  r|  tan  0]  is  the  linear  phase  distribution  in  radians  per  unit 
o 

length.  From  the  geometry  of  Figure  20, 


Figure  20.  Position  of  Slant  Range  r  for  a  Constant  Doppler 
— Frequency  Shift 


69 


t|  =  h  tan  0  +  vt 


Therefore 


d<fr 

dt 


=  r tan  6  £  [h  tan  9  +  vt] 
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which  in  terms  of  frequency  reduces  to 


w  =  w  —  tan  0 
o  c 


(7 


Equation  (78)  is  equivalent  to  the  last  term  of  Equation  (77),  for  in  both 
equations  the  frequency  shift  is  given  by  the  product  of  the  component  of 

Uq 

velocity  normal  to  the  wavefront  and  the  factor  —  . 

Hence,  it  may  again  be  concluded  that  the  frequency  of  the  signal 
received  from  an  incremental  area  centered  about  the  point  P'  is  con¬ 
stant  provided  that  r  is  normal  to  the  wavefront. 
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VI.  EXPERIMENTAL  DEMONSTRATION 


An  analytic  investigation  of  this  type  should  be  accompanied  by  an 
experimental  program  designed  to  verify  the  mathematical  description 
of  the  system  and  to  establish  the  practical  feasibility  of  the  approach. 
Generally,  this  experimental  program  can  be  of  two  types.  One  is  a  pro¬ 
gram  carried  out  in  the  laboratory,  where  modeling  or  scaling  of  the 
physical  components  is  done  to  achieve  convenient  values  of  size,  fre¬ 
quency,  etc.  The  other  type  is  an  actual  field  measurement  program 
with  a  minimum  or  breadboard  system  using  conveniently  available 
existing  facilities.  Application  of  each  of  these  avenues  of  approach  to 
an  experimental  demonstration  of  the  Transverse  Doppler  Pattern  Meas¬ 
urement  Technique  is  discussed  below,  with  emphasis  on  performing  an 
adequate  demonstration  with  a  modest  amount  of  time  and  effort. 

As  indicated  in  the  material  presented  above,  this  technique  is 
most  readily  applied  to  linear  arrays  or  rectangular  apertures  with  sep¬ 
arable  distributions.  Accordingly,  the  experimental  programs  discussed 
assume  linear  arrays  as  the  radiating  structures.  Extension  to  rectan¬ 
gular  arrays  is,  of  course,  straightforward. 

A.  A  Laboratory  Demonstration 


Since  this  program  is  concerned  with  the  measurement  of 
the  patterns  of  large  antennas  (large  both  physically  and  electrically), 
modeling  or  scaling  of  the  radiating  structure  and  the  probe  vehicle  is 
required  for  work  in  the  laboratory.  Reduction  of  both  the  physical  and 
electrical  sizes  of  these  structures  may  be  accomplished  by  scaling  in 
frequency.  If  one  somewhat  arbitrarily  chooses  100X  as  a  typical  aperture. 
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electrical  length  and  10  feet  as  a  practical  length  for  experimentation  in 
the  laboratory,  these  sizes  can  be  obtained  at  a  frequency  of  9830  Me. 
Using  Figure  3  (page  16),  for  ka  =  628.4,  a  =  2,  and  a  linear  resolution 
corresponding  to  N  =  10,  kv/B  =  17.6.  With  B  =  30  cps,  v  =  8.4  ft/sec. 
From  Figure  4,  kh  =  197,  giving  h  =  3.14  feet.  Summarizing,  the  Bys¬ 
tem  parameters  then  are: 

L  =  100X.  =  10  feet 
f  =  9830  Me 
h  =  3.14  feet 
v  =  8.4  ft/  sec 
B  =  30  cps 
N  =  10  . 

A  quite  compact  arrangement  may  be  designed  with  these  parameters. 
The  main  problem  is,  of  course,  that  any  reasonably-sized  probe  vehicle 
is  too  small  to  carry  the  required  signal  processing  system.  Therefore, 
a  means  of  transferring  the  signal  received  at  the  probe  antenna  to  a  sta¬ 
tionary  signal  processing  system  must  be  provided.  Perhaps  the  simp¬ 
lest  way  is  to  transfer  the  signal  with  a  coaxial  transmission  line  of 
sufficient  length  to  allow  the  probe  vehicle  to  move  freely  over  its  entire 
path. 


A  schematic  of  this  laboratory  system  and  a  conceptual 
drawing  of  its  layout  are  shown  in  Figures  21  and  22.  In  Figure  21  the 
linear  array  is  indicated  at  the  left  with  the  movable  probe  antenna  shown 
to  the  right  of  it.  A  portion  of  the  transmitter  signal  is  used  as  the  local 
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Figure  21.  Schematic  of  Laboratory  Experiment 


oscillator  signal  and  mixed  with  the  signal  received  by  the  probe  antenna 
to  form  a  homodyne  (or  0  IF)  receiving  system.  The  mixer  is  followed 
by  a  low-pass  filter  to  provide  the  resolution  on  the  aperture,  phase,  and 
amplitude  detector,  and  a  dual  trace  recorder. 

By  using  a  portion  of  the  transmitter  power  as  the  local 
oscillator  signed,  good  audio  frequency  stability  can  be  obtained  without 
requiring  excessive  stability  in  the  RF  signal.  Narrowband  filtering  of 
the  Doppler  component  can  then  be  accomplished  without  difficulty. 
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Figure  22.  Conceptual  Drawing  of  the  Laboratory  Experimental  System 


Figure  22  is  an  artist's  conception  of  the  experimental  setup. 
The  antenna  whose  far  field  pattern  is  to  be  measured  consists  of  a  100X. 
X-band  waveguide  linear  array.  Suspended  above  the  array  on  a  nylon 
line  is  a  small  carriage  for  supporting  the  probe  antenna.  Since  the  aper¬ 
ture  length  is  10  feet,  a  probe  path  length  of  approximately  20  to  30  feet 
will  be  required  to  adequately  simulate  actual  full  scale  measurement 
conditions.  A  high  speed  synchronous  electric  motor  with  its  shaft  speed 
reduced  in  a  gear  train  supplies  the  motive  power  for  the  carriage.  The 
motor  is  located  at  one  end  of  the  probe  carriage  path  and  draws  the  car¬ 
riage  toward  itself  with  a  nylon  cord.  A  second  nylon  cord,  attached  to 
the  opposite  end  of  the  carriage,  leads  to  a  device  which  provides  a  small 
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constant  tension  opposing  the  motion  of  the  carriage.  This  insures  that 
the  probe  velocity  will  remain  essentially  constant  over  the  major  por¬ 
tion  of  its  travel.  An  automatic  limit  switch  is  provided  to  disengage  the 
motor  as  the  extremity  of  the  carriage  travel  is  approached.  The  probe 
antenna  is  an  X-band,  half -wavelength  dipole  suspended  from  the  car¬ 
riage.  The  RF  signal  received  by  the  probe  antenna  is  fed  through  a 
length  of  flexible  miniature  shielded  coaxial  transmission  line  to  the 
signal  processing  network,  then  detected  and  recorded.  The  recorded 
trace  is  a  plot  of  received  signal  versus  time  and  corresponds  to  the 
aperture  distribution.  Both  phase  and  amplitude  detection  must  be  per¬ 
formed  in  order  to  determine  the  distribution.  The  far  field  character¬ 
istics  are  then  determined  by  integration  of  the  aperture  distribution. 

Evaluation  of  the  system’s  performance  in  measuring  vari¬ 
ous  aperture  distributions  may  be  made  by  measuring  the  aperture 
distributions  of  the  antennas  and  comparing  the  far  field  patterns  com¬ 
puted  from  the  measured  aperture  distributions  with  the  actual  patterns 
measured  by  conventional  means  in  the  far  fields  of  the  antennas.  Fab¬ 
rication  of  individual  slot  waveguide  arrays  for  a  variety  of  amplitude 
and  phase  distributions  is  required  to  obtain  the  test  arrays. 

An  alternate  arrangement  of  the  experimental  setup  is 
shown  in  Figure  23.  Here  the  probe  path  is  vertical  rather  than  hori¬ 
zontal  and  some  of  the  difficulty  in  maintaining  a  constant  distance  be¬ 
tween  the  probe  and  the  aperture  is  alleviated. 
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B. 


Field  Demonstration 


Because  of  the  availability  of  several  large  antennas  located 
throughout  the  country,  it  is  possible  that  a  full  scale  demonstration 
might  be  preferable  to  the  scaled  experiment  outlined  above.  Certainly 
in  the  final  analysis,  a  full  scale  program  must  be  carried  out;  hence  if 
the  laboratory  experiment  can  be  avoided,  much  time  and  expense  will 
also  be  avoided. 

The  apparent  advantages  of  a  full  scale  approach  then  are 

twofold: 


1.  Elimination  of  the  design  and  fabrication  of 
a  transmitting  antenna. 

2.  More  realistic  results. 

For  example,  consider  the  long  linear  array  located  near  Archer  City, 
Texas  (part  of  the  U.  S.  Navy's  space  surveillance  system);  this  array 
has  a  maximum  linear  dimension  of  one  mile.  At  the  frequency  of 
operation,  this  corresponds  to  a  length  of  the  order  of  1000  wavelengths. 
Typical  parameters  for  an  antenna  of  this  size  are: 

a  =  2 
N  =  60 

V  =  840  ft/sec 
B  =  30  cps 
h  =  1190  feet 

These  values  have  been  computed  from  the  curves  of  Figures  3  and  4. 
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Using  these  parameters,  the  experimental  procedure  differs 
from  that  of  the  laboratory  technique  in  essentially  two  ways: 

1.  A  portion  of  the  transmitter  output  can  no  longer 
be  used  as  a  local  oscillator  signal. 

2.  Variations  in  probe  path,  both  direction  and 
velocity,  must  be  considered. 

Of  the  former,  no  great  difficulty  is  encountered  since  local  oscillators 

g 

with  frequency  stability  factors  of  the  order  of  one  part  in  10  are  well 
within  the  state-of-the-art. 

Concerning  probe  path  variations,  as  has  been  stated  in 
Section  IV.  B,  errors  of  this  type  can  be  detected  and  compensated  for 
by  means  of  a  suitable  inertial  system;  however,  since  the  fly-by  time 
is  low,  it  is  felt  that  these  effects  can  be  ignored  in  the  initial  experi¬ 
mental  investigation. 

It  should  be  noted  that  all  data  processing  and  recording 
equipment  can  be  carried  by  the  mapping  vehicle,  where  the  stable  local 
oscillator  frequency  is  determined  by  the  desired  off-set  from  the  ground 
transmitter  frequency. 

4- 

Since  the  array  considered  here  is  a  discrete  element  array, 
the  aperture  distribution  is  available  by  actual  element  measurement. 
Hence  the  antenna  provides  an  accurate  experimental  standard.  Figure 
24  illustrates  the  mapping  vehicle,  the  test  array,  and  the  data  pro¬ 
cessing  system.  The  sampling  probe  is  a  single  element  which  provides 
large  angular  coverage.  The  received  signal  is  mixed  with  a  signal 
equal  in  frequency  to  the  ground-stationed  transmitter,  and  the  output 
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Figure  24.  Schematic  of  Experimental  Field  Demonstration 


of  the  mixer  is  then  passed  through  a  low  pass  filter.  The  output  of  the 
filter  then  contains  the  desired  information  concerning  the  average  phase 
and  amplitude  of  the  energy  radiating  from  the  sampling  strip. 

It  is  felt  that  data  obtained  in  an  experiment  of  this  type  would 
provide,  at  a  nominal  expense,  much  information  on  the  feasibility,  limits, 
and  accuracy  of  Transverse  Doppler  Pattern  Measurements. 
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VII.  SUMMARY  AND  CONCLUSIONS 


A  mathematical  description  of  the  Transverse  Doppler  Pattern 
Measurement  Technique  has  been  derived  on  the  basis  of  the  scalar  dif¬ 
fraction  equation.  In  the  initial  phase  the  effect  of  Doppler  frequency 
filtering  was  analyzed  and  the  results  of  this  analysis  were  applied  to 
specific  and  one-  and  two-dimensional  apertures.  The  data  computed 
from  this  application  have  established  certain  limitations  as  to  the  ver¬ 
satility  of  a  Doppler  filter  sampling  system.  Specifically,  it  has  been 
shown  that  this  technique  provides  a  sampling  resolution  in  one  dimen¬ 
sion  only  and  is  therefore  most  applicable  to  linear  arrays  and  rectan¬ 
gular  apertures  with  separable  distribution. 

In  an  attempt  to  develop  a  technique  suitable  for  more  general  ap¬ 
ertures  and  distributions,  a  second  analysis  was  performed  which  included 
the  effect  of  a  linear  array  sampling  antenna.  Because  of  the  additional 
resolving  power  of  the  directional  sampling  antenna,  it  was  possible  to 
derive  a  sampling  function  which  isolated,  within  the  boundaries  of  an 
arbitrary  two-dimensional  aperture,  a  rectangular  sampling  area  of  com¬ 
paratively  small  dimensions. 

In  an  attempt  to  extend  this  approach,  the  effect  of  multichannel 
processing  was  examined  in  which  the  radiating  aperture  was  broken  up 
into  a  series  of  incremental  areas  by  means  of  a  multichannel  Doppler 
filter  bank.  It  was  found  that  while  this  concept  is  mathematically  feas¬ 
ible,  limitations  are  again  encountered  because  of  the  directional  depen¬ 
dency  of  the  signal  received  from  each  sampling  increment.  Since  this 
directional  dependency  is  a  function  of  sampling  increment  size,  an  at¬ 
tempt  to  increase  the  sampling  system  resolution  was  made  by  utilizing 
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the  concept  of  cynthetic  aperture  processing.  Here  again  the  coherence 
of  the  quasi-point  sources  which  make  up  the  aperture  distribution  proved 
to  be  detrimental  and  it  was  shown  that  synthetic  aperture  processing  is 
not  applicable  to  the  Transverse  Doppler  Pattern  Measurement  Technique. 

It  has  been  concluded  that  the  technique  developed,  while  in  some 
respects  limited,  are  suitable  for  measuring  with  good  accuracy  a  vari¬ 
ety  of  large  one-  and  two-dimensional  arrays.  The  extent  of  the  accu¬ 
racy  is  of  course  dependent  upon  many  factors  and  has  been  discussed 
in  some  detail  in  the  Error  Analysis  section  of  the  report. 

Specifically,  it  is  felt  that  long  linear  arrays  and  two-dimensional 
arrays  with  separable  distributions  can  be  mapped  with  a  minimum  of 
equipment  and  expense.  Hence,  two  experimental  procedures  have  been 
outlined  c.i  this  basis  which  are  designed  to  prove  the  feasibility  of  the 
analytic  results. 

However,  it  is  emphasized  that  while  the  possibility  of  mapping  an 
entire  aperture  in  one  probe  transit  does  not  appear  to  be  feasible,  a 
technique  has  been  developed  which  provides  for  the  principal  plane  dis¬ 
tribution  mapping  of  any  aperture. 

In  the  case  of  a  one -dimensional  array,  the  result  is  of  course  the 
linear  distribution.  For  a  two-dimensional  array  with  separable  distri¬ 
bution,  sampling  by  means  of  a  simple  probe  produces  a  mapping  of  the 
aperture  distribution  in  the  direction  of  probe  transit  which  is  weighted 
by  the  transverse  aperture  distribution.  In  this  case,  a  second  orthog¬ 
onal  cut  is  needed  in  order  to  separate  the  mutual  weighting  effects  of 
the  distribution  functions.  If  a  linear  array  sampling  probe  is  used, 
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principal  plane  aperture  distribution  cuts  can  be  obtained  from  any  type 
of  radiating  aperture.  It  has  been  shown  that  even  the  large  phase  de¬ 
viations  associated  with  electronically-scanned  apertures  present  no 
problem. 

Hence,  it  is  apparent  that  the  Transverse  Doppler  Pattern  Meas¬ 
urement  Technique  is  feasible  and  has  great  potential  in  the  expanding 
field  of  large  aperture  antennas.  It  is  felt  that  the  next  logical  step  is 
the  development  of  an  operational  system.  As  outlined  in  Section  VI,  a 
scaled  laboratory  experiment  or  a  rudimental  field  measurement  would 
provide  much  information  concerning  the  actual  hardware  needed  to 
detect,  process,  and  record  the  sampling  system  data.  Having  accom¬ 
plished  this,  it  remains  to  refine  the  system  resulting  in  a  reliable  and 
accurate  system  for  obtaining  the  radiation  characteristics  of  large 
aperture  antennas  from  measurements  in  the  near  field. 


APPENDIX  A 


A  formulation  which  appears  to  somewhat  simplify  the  computa¬ 
tional  difficulties  of  Equation(32)  may  be  obtained  by  retaining  the  aper¬ 
ture  coordinates  in  polar  form.  Under  these  conditions,  the  equation 
corresponding  to  Equation  (8)  is 
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As  shown  in  Appendix  B,  the  t'  integral  may  be  performed  to  give 
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integral  in  u  may  be  performed  also  for 


H(m) 


sin  (w 


«  ) 
o 


(“  -  W  )  \ 


with  the  approximation 


<  <  < 


give 


|  t  -  p  sin  ip|  <  -y”! 


|  t  -  p  sin  t|»|  >  jJ 


(8 


(8 


With  the  substitution  of  Equations  (80)  and  (82),  Equation  (79)  becomes 


Sy(t)  = 


-k  exp  (jwQt) 


2vT 


r-r 


pF(p)  H 


(1) 


k  4 


.,2  2  2  n 

(h  +  p  cos  ip) 


dp 


(8 


for  |  vt  -  p  sin  >Jj| 


and 


Sy(t)  =  0  (8 

v  T 

for  |  vt  -  p  sin  ^  Equation  (83)  still  contains  a  two-dimensional 

integral,  but  it  has  somewhat  different  properties  from  Equation  (15):  The 
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sampling  function  is  not  expressed  explicitly,  but  appears  implicitly 
in  the  limits,  and  the  only  approximation  made  in  the  evaluations  is 
that 
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This  condition  is  met  very  well  for  typical  system  parameters,  such 
as  those  discussed  above. 

That  Equation  {  83)  has  basically  the  same  character  as  Equa¬ 
tion  (  15)  may  be  seen  by  considering  the  asymptotic  form  of  ^(z) 
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With  this  substitution,  the  form  of  the  two  equations  is  identical  (p  cos  ijj 

=  £)  with  the  exception  of  the  — -  —  sampling  function  and  an  amplitude 
2  2  2  U 

factor  (h  +  p  cos  ^  appears  as  the  fourth  root  rather  than  the  square 
root).  However,  as  mentioned  above,  the  sampling  dependence  is  im¬ 
plicit  in  the  limits  of  Equation  (83). 
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APPENDIX  B 


To  evaluate  the  integral 
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and  the  function 


One  way  to  solve  this  last  integral  is  to  separate  the  integral  into  two 
integrals  over  different  parts  of  the  region  of  integration. 
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give 8  the  relation 


X  =  Si.  fj  (z)  +  jY  (z)}  =  H(1)  (z) 

1  v  L  o  o  J  v  o 

where  J  (z),  Y  (z),  and  (z)  are  the  Bessel  functions  of  zero  order 

°°°  (1) 
and  argument  z  that  are  commonly  designated  by  J  ,  Y  ,  and  H 

Reintroducing  the  value  of  z  gives 
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